Centrality-dependent Direct Photon p t spectra in Au+Au Collisions at RHIC 
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We calculate the centrality-dependence of transverse momentum (p t ) spectra for direct photons 
in Au+Au collisions at the RHIC energy, based on a realistic data-constrained (3+1) dimensional 
hydrodynamic description of the expanding hot and dense matter, a reasonable treatment of the 
propagation of partons and their energy loss in the fluid, and a systematic study of the main sources 
of direct photons. The resultant pt spectra agree with recent PHENIX data in a broad pt range. The 
competition among the different direct photon sources is investigated at various centralities. Parton 
energy loss in the plasma is considered for photons from fragmentation and jet photon conversion, 
which causes about 40% decrease in the total contribution. In the high p t region, the observed Raa 
of photons is centrality independent at the accuracy of 5% based on a realistic treatment of energy 
loss. We also link the different behavior of Raa for central and peripheral collisions, in the low p t 
region, to the fact that the plasma in central collisions is hotter compared to peripheral ones. 

PACS numbers: 



I. INTRODUCTION 

The formation and observation of a quark-gluon 
plasma (QGP) in heavy ion collisions are important goals 
of modern nuclear physics [l|, 0| ■ Suppression of high p t 
hadron yields [3( is one of the most important features 
observed at the Relativistic Heavy Ion Collider (RHIC). 
Theoretically this is attributed to the interaction between 
jets (hard partons) and the bulk matter [J, [f|, [y, 0] . Ex- 
perimentally, absence of the suppression in <i+Au colli- 
sions [8| reveals that the suppression results from a final- 
state effect and, in turn, that the hot and dense matter is 
created in Au+Au collisions. The amount of suppression 
depends significantly on the centrality of the collision Q , 
which implies various sizes of hot dense matter are formed 
in heavy ion collisions at various centralities. This offers 
us an excellent opportunity to study the interaction of 
partons inside the system and, consequently, properties 
of the matter under extreme conditions. 

Hadron production in heavy ion collisions involves bulk 
hadronization of the thermal partons at low pt, the frag- 
mentation of quenched hard partons at high p t and the 
hadronization contributed from both thermal and hard 
partons at intermediate pt- However, it is quite difficult 
to systematically describe all these hadronization pro- 
cesses since some of them are beyond the perturbative 
treatment and usually contain many parameters without 
full understanding. Low p t hadrons also strongly inter- 
act with each other after hadronization and cannot carry 
direct information from inside the hot matter. Under 
this situation, a systematic study of direct photons in a 
wide range of transverse momentum and centrality can 
serve as a guide to understand the whole reaction pro- 
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cesses of heavy ion collisions, since we do not need to 
treat hadronization itself nor interaction between pro- 
duced direct photons and the bulk matter thanks to the 
large mean free path of direct photons compared to the 
typical size of the system in heavy ion collisions. Compe- 
tition among different sources at various centralities may 
be also useful in understanding of production mechanism 
of direct photons. 

In this paper, we first study the role of jet quenching 
on the centrality dependence of direct photon production. 
For this purpose, a reliable treatment of hard parton en- 
ergy loss is needed. This is formulated via the BDMPS 
framework [7J and tested on pion suppression at various 
centralities. Since neutral pions and other mesons are 
significantly suppressed in central Au+Au collisions, and 
since the suppression has an evident centrality depen- 
dence, the following question arises naturally: What is 
the role of hard parton energy loss on direct photon pro- 
duction? Main purpose of this paper is to answer this 
question. 

We also investigate the interplay among the various 
sources of direct photons. Similar to hadron suppres- 
sion, photons from parton fragmentation are expected to 
offer information on the interaction between hard par- 
tons and the bulk via jet quenching. Thermal photons 
and photons from parton-bulk interactions are penetrat- 
ing probes of the hot matter, respectively, through in- 
teraction of partons inside bulk matter and interaction 
between primary partons and the bulk matter. It is in- 
teresting to see whether one reproduces the observed pho- 
ton spectra by considering all photon sources simultane- 
ously and consistently at different collision centralities. 
By identifying the dominant sources of direct photons 
at given values of pt, we will be able to discuss in which 
way the different p t -regions of the photon spectra provide 
information about the different production processes. 

We need a realistic description of the hot and dense 



matter to investigate the effect of bulk matter on photon 
emission. This is achieved by using three dimensional 
(3D) hydrodynamic simulations of bulk matter [lfj, [jj| 
which have already been tested against a vast body of 
low p t hadron data at RHIC. 

To perform a systematic study of direct photon pro- 
duction (from sources other than neutral meson decays) 
in relativistic heavy ion collisions, we shortly review the 
possible sources in the following. 

Primordial NN scattering. The direct photon pro- 
duction via Compton scattering and quark-antiquark an- 
nihilation can be calculated in perturbation theory us- 
ing the conventional parton distribution functions and 
the factorization hypothesis. In principle one should 
consider at this stage also higher order contributions, 
like bremsstrahlung of photons accompanying for exam- 
ple two-jet production in hard parton-parton scattering. 
However, we consider this component as a part of the 
so-called jet fragmentation (or bremsstrahlung) contri- 
bution to be affected by the thermalizcd matter, which 
we will discuss separately. 

Thermal photons. In high energy nuclear collisions, 
the density of secondary partons is so high that the 
quarks and gluons rescatter and eventually thermalizc 
to form a bubble of hot QGP. The plasma expands, cools 
down, and goes through a phase transition to hadronic 
gas (HG) phase. Thermal photons can be produced dur- 
ing the whole history of the evolution of hot matter from 
the QGP phase to the HG phase through the mixed phase 
due to collisions of or radiations from thermalizcd par- 
ticles. Yields of photons from a thermal source are ex- 
ponentially damped so that contribution to very high p t 
region is negligible. However, contribution to low p t is 
expected to be dominant in central collisions in which 
the size and the temperature of a hot and dense matter 
are large enough. 

Jet-photon conversion. When hard partons pass 
through thermalizcd matter, they may interact. Col- 
lisions between jets and deconfined partons via quark- 
antiquark annihilation and quark-gluon Compton scat- 
tering can produce direct photons. This is often called 
as jet-photon conversion. 

Jet fragmentation. Photon production also occurs 
as a higher order effect in purely partonic initial hard 
scatterings: at any stages of the evolution of a jet (final 
state parton emission) , there is a possibility of emitting 
photons. Existence of a QGP again affects the results 
of fragmented photons since energetic partons lose their 
energy prior to fragmentation. In this work, we assume 
fragmentation of partons only outside the plasma, which 
is similar to high p t hadron production from jet fragmen- 
tation. 

There are possible contributions to photon produc- 
tion which are not included in the present study: The 
medium-induced radiation is supposed to be a higher or- 
der contribution. At the RHIC energy this contributes 
much less than jet photon conversion at low and inter- 
mediate p t and much less than fragmentation at high 



Pt [12J |. So the contribution from medium-induced ra- 
diation is ignored this paper. In the time interval be- 
tween the primordial collisions at r = and the thermal- 
ization of the hot matter at To, the interaction between 
non-equilibrated soft partons and hard partons may also 
produce direct photons. We neglect the contribution in 
the prcequilibrium stage since the time interval is much 
shorter than the life time of the equilibrated matter (~ 20 
fm/c) . 

The paper is organized as follows: In Sec. HH we first 
give a brief review on the space-time evolution of the hot 
matter created in Au+Au collisions at different ccntrali- 
ties based on a (3+l)-dimensional ideal hydrodynamical 
calculation. In Sec. IIII1 we discuss parton energy loss 
in the QGP. We investigate neutral pion production in 
the high p t region in order to fix the parameters of the 
energy loss scheme. We discuss sequently the contribu- 
tions from various sources to direct photon p t spectra in 
Scc. lIVI We show our results and compare them with re- 
cent experimental data in Sec. [V] Section IVll is devoted 
to conclusion of the present study. 



II. SPACE-TIME EVOLUTION OF THE HOT 
AND DENSE MATTER 

Several sources of direct photon production in heavy 
ion collisions depend on the bulk dynamics of hot and/or 
dense matter and the matter along trajectories of en- 
ergetic partons. So a realistic description of reaction 
dynamics is indispensable for the quantitative analysis 
of photon production. In our calculation, fully three- 
dimensional (3D) ideal hydrodynamics [l(J, LUj is em- 
ployed to describe the space-time evolution of the hot 
and dense matter created in Au+Au collisions at RHIC 
energy at various centralities. We solve the equations of 
energy-momentum conservation 



duT"" 







(1) 



in full 3D space (r,x,y,r]) under the assumption that 
the local thermal equilibrium is reached (maintained) at 
(after) an initial time To =0.6 fm/c. Here r and rj are the 
proper time and the space-time rapidity, respectively, x 
and y are transverse coordinates. In the transverse plane, 
the centers of two colliding nuclei are located at (x, y) = 
(6/2,0) and (—6/2,0) before the collision at an impact 
parameter b. Ideal hydrodynamics is characterized by 
the energy- momentum tensor, 



T» v = (e + P)uV - Pg" 



(2) 



where e, P, and w^ are energy density, pressure, and lo- 
cal four velocity, respectively. We neglect the finite net- 
baryon density which is small near the mid-rapidity at 
RHIC. For the high temperature (T > T c = 170 MeV) 
QGP phase we use the equation of state (EOS) of mass- 
less non-interacting parton gas (u, d, s quarks and glu- 



Table I: Initial temperature at the plasma center at initial 
time tq = 0.6 fm/c for various centralities. 



Centrality(%) 


0-10 


10-20 


20-30 


30-40 


40-50 


50-60 


T (MeV) 


370 


357 


341 


327 


301 


272 



ons) with a bag pressure B: 



1 



p=-(e-4B). 



(3) 



The bag constant is tuned to be S" = 247.19MeV to 
match pressure of the QGP phase to that of a hadron 
resonance gas at critical temperature T c = 170 MeV. 
A hadron resonance gas model at T < T c includes all 
hadrons up to the mass of the A(1232) resonance. Our 
hadron resonance gas EOS implements chemical freeze- 
out at T c h = T c = 170 MeV, as observed in collisions 
at RHIC [13| ■ This is achieved by introducing appro- 
priate temperature-dependent chemical potentials fJ-i(T) 
for all hadronic species i in a way that their numbers N 
including all decay contributions from higher-lying reso- 
nances, Ni = N + y] R bR^ixNjj, are conserved during 
the evolution [ll|, [lj, jjj, |lj, M Ell • Hcrc N i is the av- 
erage multiplicity of the i-th hadron species, and &_r^a" 
is the effective branching ratio (a product of branching 
ratio and degeneracy) of a decay process R — > i + X. In 
this partial chemical equilibrium (PCE) model [ll| only 
strongly interacting resonances with large decay widths 
(whose decays do not alter N) remain chemically equili- 
brated below T c h. It should be noted that the hadronic 
chemical composition described by hydrodynamics using 
the PCE model is roughly consistent with that of the 
hadronic cascade models [lj|, as long as the latter are 
initialized at T sw = 169 MeV with thermal and chemical 
equilibrium distributions. 

We assume that, at tq = 0.6 fm/c, the initial entropy 
distributions is proportional to a linear combination of 
the number density of participants (85%) and binary 
collisions (15%) [19|. Centrality dependence of charged 
particle multiplicity observed by PHOBOS [2(| has been 
well reproduced by full 3D hydrodynamics simulations 
with the above setups [lj|. In the following calculations, 
hydrodynamic outputs at representative impact parame- 
ters b = 3.2, 5.5, 7.2, 8.5, 9.7, and 10.8 fm are chosen for 
0-10%, 10-20%, ■ ■ • , 50-60% centrality, respectively. 

So far, the space-time evolution of the QGP fluid ob- 
tained as above has been also exploited for a quantita- 
tive study of hard and rare probes such as azimuthal 
jet anisotropy, nuclear modification factor of identified 
hadrons, disappearance of back-to-back jet correlation, 
and J/tp suppression [211 ]. 

In Table HI initial temperatures at the plasma cen- 
ter, T = T(t , 0, 0, 0), are shown for various centralities. 
These temperature values will be important to interpret 
the centrality dependence of the slope of pt spectra from 
thermal radiation, which will be discussed later. Fig- 
ure Q] shows the time evolution of energy density at the 
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Figure 1: (Color Online) Time evolution of energy density at 
the center (x, y, rj) — (0, 0, 0) for various centrality. Each line 
from top to bottom corresponds to 0-10%, 10-20%, ..., and 
50-60% centrality, respectively. 



center of fluids (x, y, rj) = (0, 0, 0) for various centralities. 
Clearly for any given proper time r, the more central col- 
lisions one obtains higher energy densities at the plasma 
center. 

For convenience of the following calculations, we intro- 
duce /qgp( t j x i V, v) as the fraction of the QGP phase in 
a fluid element. It is obvious that /qgp = 1 (0) in the 
QGP (hadronic) phase. In the mixed phase, the fraction 
of the QGP is calculated via 

,/qgp e QG p + (1 - ./QGp)ehad = e(r, x, y, 77) 

with eQGP and ehad being the energy densities of the 
QGP phase and the hadron phase at T — T c , respectively. 



III. PARTON ENERGY LOSS IN A PLASMA 

Energy loss of hard partons in a plasma affects both jet 
photon conversion and jet fragmentation. The momen- 
tum distribution of jets (energetic gluons or quarks with 
different flavors) from primordial nucleus-nucleus scat- 
tering is calculated as |22| 



dN AB^jet 

dyd 2 pt 



abed 



= KT AB {b) J2 / dx a dx b G a/A (x a ,M 2 ) 



8 da 



xG b/B (x b , M 2 )--,{ab -► cd)6(s + t + u) (4) 

7T dt 

where Tab (b) is the nuclear overlapping function at an 
impact parameter b for each centrality, G a / A (x a , M 2 ) and 
G b /g(xb, M 2 ) are parton distribution functions in nuclei 
A and B. We take MRST 2001 LO parton distributions 
in proton [23j . The elementary cross sections for ab — > cd 
can be found in Ref. [221 . We set the factorization scale 



M and renormalization scale Q to be M = Q = pt- 
K = 2 is chosen to take into account higher order con- 
tributions. These parameters are chosen as to reproduce 
high pt pion data in pp collisions at RHIC, which will 
be discussed later. The above formula for p t spectra was 
extensively tested in pp (pp) collisions in an energy range 
from y/s = 27 A GcV to 630 GcV. Nuclear shadowing 
effect and EMC effect are taken into account through 
EKS98 scale dependent nuclear ratios R^ KS (x, A) [24| . 
Isospin of a nucleus with mass A, neutron number N, 
and proton number Z is corrected as follows: 



G a /A(x) 



N 'n , s Z r, / x 

-jG a/N (x) + —G a/P (x) 



R 



EKS 



(x,A). (5) 



The isospin mixture and nuclear shadowing eventually 
cause a decrease of nuclear modification at high p t region, 
which will be shown in Sec. [Vj 

We assume that all jets are produced at r = 1/Q w 
with the phase space distribution 



. . _ . dN ( b 
fo(p, r) ex -^~T A \x--,y)T B [x 



y S(z) (6) 



where r = (x, y, z) is the coordinate of a jet, b is the 
impact parameter, and Ta and Tg are thickness functions 
of nuclei A and B. The <5-function reflects the highly 
Lorentz-contracted colliding nuclei A and B. The phase 
space distribution of hard partons is normalized as 



JM^r = (2tt) 3 |^ 



(7) 



Energetic partons can suffer interactions with the fluid 
and lose their energies. We employ the BDMPS for- 
mula [y] to calculate parton energy loss in a plasma cre- 
ated in heavy ion collisions. For a parton of type i = q, g 
with initial momentum p formed at r , the whole path 
length of a parton traversing the expanding QGP (in- 
cluding the mixed phase) is 



L(p ,r ) 



dr9(f QGP (T,x(r))). 



(8) 



Here x(t) is a trajectory of a parton, /qgp(t, x(t)) is the 
fraction of the QGP phase at a position (t,x(t)), and 9 
is a step function, which gives 6(j qgp) equal unity in 
the QGP and the mixed phases and zero in the hadron 
phase. 

The total energy loss along this path is calculated as 

AE(i,p ,f ) = D / dTe(i,T,x(r))e(f QGP (T,x(T))). 

J To 

(9) 
Here D is an adjustable parameter, e(i, r, x(t)) is the 
energy loss per unit distance for a parton i at a position 
(t,x(t)), given as [|| 

e(i, t, x(r)) = a s \/fi 2 E*/\i. 



Here, the temperature-dependent running coupling con- 
stant - assuming a similar formula as the lowest order one 
in perturbation theory - can be obtained by fitting the 
numerical results from lattice quantum chromodynamics 
(QCD) simulations [25( as 



«. CO 



(J7T 



(33 - 27V/) ln(8T/T c ) 



(10) 



The Debye screening mass is given as /i = gT, with 
g 2 /47r = a s (T). The energy of a hard parton in the 
local rest frame is E* = p^u^ where p^ is the four mo- 
mentum of the hard parton in the laboratory frame and 
Ufj, is a local fluid velocity. All hard partons are treated 
as on-shell massless particles. The mean free path, A;, of 
a parton i, is given as 

K = VgqPqfQGP + VggPgfQGP, (H) 

K = VqqPqfQGP + OqgPgfo,GV, (12) 

with cross sections <Zj = Cia s n/T 2 [2y|. The color fac- 
tors 2d are 4/9, 1, and 9/4 for qq, qg, and gg scattering, 
respectively. The parton densities p q and p g are obtained 
from the EOS of the massless relativistic gas. The frac- 
tion of the QGP phase /qgp is considered in the mixed 
phase to ensure a smooth transition from the QGP phase 
to the HG phase. Note that the above quantities, i.e., 
temperature T, fluid velocity « M , parton densities pi and, 
in turn, mean free path Ai, depend on the location of the 
parton x(t) and can be obtained from full 3D hydrody- 
namics simulations discussed in the previous section. 

Various sizes of the plasma arc formed in heavy ion 
collisions at different centralitics. We use the common 
energy loss formula Eq. ^ for all of these media. The 
main purpose in the present paper is a systematic study 
of direct photon production rather than a detailed anal- 
ysis of parton energy loss. So we admit ourselves to in- 
troduce an adjustable parameter D in Eq. ([9]) to fit si- 
multaneously neutral 7r-meson data in Au+Au collisions 
at different ccntralities 0. 

We first discuss pion production in proton-proton col- 
lisions. We calculate neutral 7r-mcson production assum- 
ing pQCD factorization, Eq. (TJJ, 



dyd 2 p t 



E 



dz. 



dNPP- 



1 



dyd 2 p\ 



D° 0/c (z c ,Q 2 ), (13) 



where D® ,(z c ,Q 2 ) is ir° fragmentation functions pa- 
rameterized by Kniehl et al. [27|- In Fig. [21 p t spectra 
for neutral pions in pp collisions at yfs = 200 GeV cal- 
culated with M = Q = 2pt, pt, and Pt/2 are compared 
to PHENIX data [28[ . In the high p t region where the 
pQCD is expected to work, we reasonably reproduce the 
experimental data with the above setup with K = 2 and 
M = Q = pt- We use the p t spectrum as a reference 
spectrum in the following calculations. 

The effect of parton energy loss is taken into account 
through the medium modified fragmentation function 



% * 

o 10 

X> 

S -4 

>, -6 

5 10 
e 

■° -8 

10 

-10 

10 
( 




p+p -> 7i° Vs=200GeV 

t — Q=P t /2 
% — Q=p, 
^x ■■ Q=2 Pt 

'^i, • PHENIX 


) 


5 10 15 20 

p t (GeV/c) 



Figure 2: (Color Online) Neutral pion production in pp col- 
lisions at y/s - 200 GeV is compared to PHENIX data [2J |. 
Three lines from top to bottom correspond to Q = pt/2, pt, 
and 2p t respectively. 



D v o/ c (z c , Q 2 , AE C ) which describes suppression of neu- 
tral pion yields as 



dN 



AB 



dyd 2 p t 



with m 



ZJ / dZc ru.^c ■^D v0/c (z c ,Q 2 ,AE c ). 



dyd 2 Pt z 2 



(14) 



D„o /c (z c ,Q 2 ,AE c ) 



L z' 



Dl 0/c {z' c ^) + -fDl 0/g (z' g ^) 



+e~ D° B/c (z c ,Q 2 ) 



(15) 



Figure [3] shows the nuclear modification factors for neu- 
tral pions in Au+Au collision at ^/snn = 200 GeV for 
different centralities. Solid lines arc results with an en- 
ergy loss parameter D = 1.5 and plots are the PHENIX 
data j9|. With a common value of the parameter D = 1.5, 
we can reasonably reproduce the it yields in the high 
Pt region at all centralities simultaneously. It should be 
noted that, in PHENIX data [§], there are ~ 10% nor- 
malization error and ~ 7—16% errors (depending on cen- 
trality) due to N co \\, which are omitted in Fig.[3J In the 
region p t < 5GeV/c, our results undershoot the exper- 
imental data due to absence of neutral pion production 
from bulk components in this calculation. Notice that 
low pt pion data have already been described well [l9( by 
using hydrodynamic simulations employed in the present 
study. In the following photon calculations, we always 
use the BDMPS energy loss formula ((9]) with D = 1.5. 



IV. THE DIFFERENT SOURCES OF DIRECT 
PHOTON PRODUCTION 

Leading order contribution. Similar to Eq. ([?]), the 
leading order contribution to direct photon production 
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Figure 3: (Color Online) Nuclear modification factors of tv°. 
Solid lines are calculated with the BDMPS energy loss formula 
with amplified parameter D = 1.5 (See text for details). Plots 
are PHENIX data @|. 



in nucleus-nucleus collisions reads 

dN AB^~, 



dyd 2 p t 



TAB{b)^2 / dx a dx b G a/A (x a ,M 2 ) 
ab 
da 



xG b/B {x b , M 2 )--,(ab -» 7 + X)S(s + t + u) (16) 

7T dt 

where the elementary processes ab — > 7 + X arc Compton 
scattering qg — > jq and annihilation qq — > 97. 

Fragmentation contribution. Higher order contribu- 
tions in pp collisions are due to jet fragmentation. We 



can calculate them as 



dNg* 

dyd 2 p t 



E 



dN pp ^ c 1 
^^^-^/c^O 8 ). (17) 



with photon fragmentation functions D° , c (z, Q 2 ) being 
the probability for obtaining a photon from a parton c 
which carries a fraction z of the parton's momentum. So 
pjr = pt/z c is the transverse momentum carried by the 
parton c before fragmentation and d 3 p/E = z 2 d 3 p c /E c . 
The effective fragmentation functions for obtaining pho- 
tons from partons can be calculated perturbatively. We 
use the parameterized solutions by Owens [23 ]. 

In case of heavy ion collisions, parton energy loss in a 
plasma should be taken into account. This can be done 
via modified fragmentation functions [291 ] 



D l/c (z c ,Q 2 ,AE c ) 



Z ^D°(z> c ,Q 2 ) 



Zc 



L z' 



+e-* D° /c (z c ,Q 2 ), 



with z' c = p t /{Pt - AE C ) and z' g = (L/X c )p t /AE C being 
the rescaled momentum fractions carried by the parton c 
and the emitted gluons before fragmentation. A c is mean 
free path of the parton c in the plasma, L is the path 
length of each parton traversing the plasma defined in 
Eq. |5]). Thus, in heavy ion collisions, contributions from 
fragmentation become 



dyd 2 p t 



E 



dz 



dN 



AB- 



1 



dyd 2 p c t zl 



D l/c (z c ,Q 2 ,AE c ). 



(19) 

The above formula counts only fragmented photons 
outside the plasma. In principle, when fragmentation 
into photons happens inside the plasma photon can es- 
cape the plasma due to the long mean free path. However 
it is not evident when and where fragmentation happens. 
Thermal production. The emission rate of photons is 
r = Ed 3 R/d 3 p, where R is the number of photons emit- 
ted from a medium per unit space-time volume with tem- 
perature T. Total yields of thermal photons can be ob- 
tained by summing the emission rate over the space-time 
volume as 



dN 



thermal 



dyd 2 p t 



= f d i xT(E*,T) 



(20) 



P^Uf, 



being the 



with d A x = rdrdxdydr] and E* 
photon energy in the local rest frame. Here, p 
(pt cosh y, p t cos 4>, Pt sin </>, p t sinh y) is the photon's four 
momentum in the laboratory frame and w M is a local fluid 
velocity. In our calculations, the thermal photon emis- 
sion rate covers both contributions from the QGP phase 
and the hadronic phase 

T(E*,T) - / QGP rQ GP ^(£*,T) 

+ (l-/ QGP )r HG ^(iT,T), (21) 



where /qgp and T are the fraction of the QGP phase 
and temperature of the fluid, respectively, both being 
obtained in the hydrodynamic simulations. In the above 
formula, we calculate thermal photon production above 
the thermal freeze-out temperature T t h = 100 McV. The 
photon emission rate from 2 — ► 2 processes between ther- 
mal partons, i.e., the QCD Compton process qg — ► ^q 
and annihilation qq — > 57, was first calculated with the 
hard thermal loop resummation technique [30j,[3l|. Later, 
Landau-Pomeranchuk-Migdal (LPM) interference effect 
for emitted photons turned out to be important [32| , lead- 
ing to 



N, 



d \ 2 aas „ 2 1 



rQOP-n ( ^T)-:5:(2)-pr»— 

i— 1 

x[\n(^-)+hn(2x) + C 22 (x) 



-Cbrcms(x)+a„n(x)], (22) 



(18) with x = E*/T and 



C 22 (x) = ^^- - 0.3615 + l.Ole" 1 - 35 *, (23) 



^brcmsi^j T" (-'ann^J 



0.633a; 



-1.5 



In (12.28 +l/x) 



0.154a; 



(l + a;/16.27) - 5 



In the calculation we take Nf = 3 and a temperature 
dependent running coupling as in Eq. (jTUJ) - 

Thermal photon emission in the hadronic phase re- 
sults from interactions such as tttt — + fry, irp — > ir~f and 
p — > 7T7T7, etc. Interactions of mesons or baryons with 
strangeness can also produce photons, but the contribu- 
tions are relatively small due to the phase space suppres- 
sion resulting from their heavier masses. In our work, 
photon emission rate from the hadronic phase is based on 
massive Yang-Mills (MYM) calculation j33[, where pho- 
ton production from mesons with strangeness has been 
included as well as the axial meson ai as an exchanging 
particle for non-strange initial states. Hadrons are com- 
posite objects, so form factors are considered to simulate 
finite hadronic size effects [34j ]. 

Jet photon conversion with jet energy loss. When hard 
partons propagate in a plasma, they also collide with 
thermal partons and produce direct photons via Comp- 
ton process and the quark-antiquark annihilation. We 
call this process jet-photon conversion, since it is a con- 
version of a jet into a photon with almost the same mo- 
mentum as the one of originating jet parton. Contribu- 
tion from the jet-photon conversion is calculated by inte- 
gration of conversion rate over the space-time evolution 
of the hot and dense matter in the QGP phase 



dN^ c 
dyd 2 p t 



T^(E*,T)f QGP (x,y, V ,T)d 4 x. (24) 



The photon production rate by annihilation and Comp- 
ton scattering of hard partons in the medium can be ap- 
proximated as [33, [3|| 



rJ pc ( £* r) = ^£ e 2 /g(p>)T2 



AEZT 
In— 2--C 



(25) 
where E* is the photon energy in the local rest frame, 
C = 2.323, m 2 h = g 2 T 2 /6, and the strong coupling 
a s = g 2 /iTT being temperature dependent as in Eq. (fT0|) . 
a is the electromagnetic couplings, e q and f q (p,x) are 
the electric charge and the phase-space density of a hard 
parton of flavor q. The phase space distribution of hard 
partons at r is obtained by considering parton energy 
loss as 

f(f,x) = f(p,r,r) 

= / d 3 p a f a (po,r- vt)5(p -p- vAE) 

where /o (p, r) is the phase space distribution at r = 
described in Eqs. ([6|) and ([7]). The (5- function expres- 
sion reflects the energy loss of a parton moving along a 
straight line trajectory with v = p/E = pq/Eq. AE is 
the energy loss from To to r and calculated similar to 
Eq. ([9]) but replacing the upper limit of integral 00 with 

T. 



V. RESULTS AND DISCUSSION 

In Fig. 21 the calculated p t spectra of direct photons in 
Au+Au collisions at ^/s^n = 200 GcV at centralityO- 
20% and 20-40% are compared to PHENIX data [37|,|38[. 
Here we sum over all contributions discussed in the previ- 
ous section. The theoretical results for 0-20% centrality 
are obtained as a mixture of the calculations for 0-10% 
and 10-20% centrality with a weight of 50% each; a cor- 
responding procedure applies for the 20-40% centrality 
results. The PHENIX data are reproduced within our 
multi-component model remarkably well. 

In Fig. [5] we show a detailed comparison of the calcu- 
lated pt spectra of direct photons with PHENIX data [39| 
for the centralities 0-10%, 10-20%, ..., 50-60%. Again, 
our results agree with data very well in a broad range of 
Pt and centrality. 

Since all the curves are almost parallel to each other, 
one gets more insight by using the nuclear modifica- 
tion factor Raa, obtained by dividing a pt spectrum in 
nucleus-nucleus collisions by the -/V co ii-scaled p t spectrum 
in pp collisions. In Fig. [6[ we show the invariant dif- 
ferential cross section of direct photons in pp collisions. 
The calculation includes the leading order contribution 
plus fragmentation contribution, using a scale Q = p t . 
PHENIX data are shown as open circles [33] and filled 
circles [381 ] . In high p t regions, our result agrees with the 
data reasonably well: So we use it to calculate nuclear 
modification in Fig. [7] and Fig. [TU] It also provides a 



baseline calculation with the LO contribution and frag- 
mentation contribution in Au+Au collisions. Whereas, 
in low pt regions where pQCD is not expected to work, 
our result undershoots the data slightly although the er- 
ror bars are large in data. The dashed line is a fit to the 
measured differential cross section of direct photons in 
pp collisions at the RHIC energy 



da 



dycPpt 



= 0.01834 1 



1.432 



-3.27 



ib/GeV 2 



which is employed to calculate the nuclear modification 
factor from thermal contribution in Fig. [TT] (a). 

Figure [7] shows how the nuclear modification factor 
for direct photons, Raa-, depends on centrality and on 
energy loss. Data for 0-10% centrality are taken from 
Refs. [39| and (4fJ. Figure EJa) shows centrality depen- 
dence of Raa compared to the PHENIX data. The three 
curves are respectively 0-10% (dotted line), 20-30% (solid 
line), and 40-50% (dash-dotted line). Raa has a weak 
centrality dependence at high pt region. This result is 
consistent with the observed phenomenon (39| that the 
Pt-integrated (for p t > 6 GeV/c) Raa of direct photons is 
almost independent of collision centrality. Does this im- 
ply a very weak effect from jet quenching? Figure [7f b) 
answers this question (here for the most central colli- 
sions): Comparing calculations with (dotted line) and 
without energy loss (solid line), one finds a difference of 
up to 40%. So the effect of parton energy loss is quite 
visible in the pt range between 4 GeV/c and more than 
20 GeV/c. If we would do the Raa calculations without 
energy loss, the difference between central and semipe- 
ripheral collisions would be about 20%, wheras the com- 
plete calculation gives the same result for all centralities, 
within 5%. We have to admit that we talk about small 
effects, requiring experimental data with relative errors 
of less than 5to observe the effects. 

To understand the above results, we look more closely 
into the different contributions. Parton energy loss in 
the plasma suppresses the fragmentation contributions 
and jet-photon conversion. So we study the ratios of the 
contribution with energy loss to the one without energy 
loss, as shown in Fig. [H] ((a) for fragmentation and (b) 
for jet-photon conversion). Energy loss in the plasma de- 
pends on the path length of the hard parton inside the 
plasma, which turns out to depend on the collision cen- 
trality. We do see a similar centrality dependence of the 
suppression for ir° (jet quenching effect) in fragmentation 
contributions and jet-photon conversion. 

To understand how these energy loss features affect 
the total contribution, we investigate the competition 
from different sources in Fig. [9[ for the three centrali- 
ties 0-10%, 20-30%, and 40-50%. The leading order con- 
tribution (LO) from primordial elementary scatterings 
is plotted as dotted lines, thermal contribution as dash- 
dotted lines, fragmentation contribution as dashed lines, 
and jet photon conversion as solid lines. The latter two 
are calculated with parton energy loss in the plasma (left 
plots) and without (right). For all centralities, thermal 
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Figure 4: (Color Online) Direct photon production in Au+Au collisions at centrality 0-20% and 20-40%. PHENIX data are 
shown as open circles [37| and filled circles [381 ] - 
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Figure 5: (Color Online) Direct photon production in Au+Au 
collisions at y/SNN =200 GeV for different centralities (0-1 0%, 
10-20%, ..., 50-60%). Data are obtained by PHENIX [3J|. 



photons dominate at low transverse momenta and they 
are insignificant in the high p t region. The leading or- 
der contribution from primordial elementary scatterings 
dominates in the high p t region. This contribution is 
independent of bulk volume. 

Let us first discuss the central collisions. Here fragmen- 
tation and conversion are an order of magnitude smaller 
that the LO contribution. But from Fig. [51 we know that 
is due to the strong energy loss effect. Without this en- 
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Figure 6: (Color Online) Invariant differential cross section 
of direct photons in pp collisions at ^/s — 200 GeV. Lines are 
leading order calculations plus fragmentation contributions, 
with a scale Q — pt. PHENIX data are shown as filled circles 
[H and open circles ]g$. Dashed line is a fit to PHENIX 
data. See text for details. 



ergy loss, these two contributions would be much bigger, 
and this is why we find a 40% difference between the total 
contribution with and without energy loss. 

For peripheral collisions, without energy loss, the rel- 
ative contribution from conversion is smaller compared 
to central scatterings, since the the plasma regions is 
smaller. But then also the suppression from energy loss is 
smaller for the peripheral compared to central collisions. 
So at the end, including a proper energy loss treatment, 
for both central and peripheral collisions, conversion is 
roughly an order of magnitude smaller than the LO con- 
tribution, see Fig. [UJc). The relative contribution from 
fragmentation, without energy loss, is comparable in cen- 
tral and peripheral collisions, however, the energy loss is 
smaller in the latter ones. So fragmentation is somewhat 
more important in peripheral compared to central, when 
energy loss is considered, as can be also seen from Fig.lTUl 
where the contribution to Raa from fragmentation and 
conversion is shown for the different centralities: Convcr- 
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Figure 7: (Color Online) The nuclear modification factor of direct photons in Au+Au collisions Raa- Data for 0-10% centrality 
are from Ref. ^ and Ref. [I}. (a): Raa at centrality 0-10% (dotted line), 20-30% (solid line), and 40-50% (dash-dotted line) 
respectively, (b): Raa at 0-10% centrality with energy loss (dotted line) and without energy loss (solid line). 
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Figure 8: (Color Online) The ratio of the contribution with energy loss to the one without, in fragmentation (a) and jet-photon 
conversion (b). 



sion contributes roughly 4%, for all ccntralitics, fragmen- 
tation between 5% (central) and 10% (peripheral). 

At the end, Raa is nearly centrality independent as 
shown in Fig. [Tfa) , but a realistic (and strong) partonic 
energy loss is needed in order to get this scaling behavior. 

In the low p t region, contribution from thermal ra- 
diation is of significant importance. We check the 
ccntrality-dcpcndence of the thermal contribution to 
Raa in Fig. [Tl](a). At very low p t , i.e., p t < 1 GeV/c, 
the thermal contributions to Raa at different ccntralitics 
coincide with each other. However the slope of i?^ rmal 
changes and the dominant pt region of thermal photons 
becomes smaller when one moves from central to periph- 
eral collisions. This reflects the fact that the temperature 
in the core region depends on the collision centrality as 
shown in Table [I] in Sec. HI] 

So from the thermal source, the Raa for central colli- 
sions exceed more and more the Raa for peripheral col- 
lisions, which translates into a slight overshooting of the 
central total Raa compared to the peripheral one, as seen 
in Fig. in the region pt < 4 GeV/c. 

The fractions of thermal contribution as a function of 



p t from different phases are shown in Fig. [TT] (b) . Par- 
tial chemical equilibrium in the hadronic phase is used 
in this hydrodynamic simulation to keep the number of 
hadrons fixed below T c h. If we ignore the contribution 
from particle decays and use a full chemical equilibrium 
(FCE), the photon emission rate from hadronic phase 
[33j can be used in this case. In case of the PCE, the 
chemical potential /x, for all hadronic species i will mod- 
ify the photon emission rate from hadronic gas, roughly 
estimated by a factor of exp[(/ii + fi2)/T] for a subpro- 
cess of 1 + 2 ^ 3 + 7 according to kinetic theory with 
a Maxwell-Boltzmann statistics for all particles. This 
factor finally increases contribution from the HG with a 
factor about 2. Nevertheless, in the total thermal con- 
tribution, PCE or FCE consideration does not make a 
visible difference. For all centralities from 0-10% to 50- 
60%, the QGP phase emits most of direct photons above 
Pt ~1 GeV/c. Although the volume of hadronic phase is 
much bigger than the one of the QGP phase due to expan- 
sion, the photon emission rate from the hadronic phase 
decreases even faster with temperature. The competition 
between volume and emission rate results in the biggest 
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Figure 9: (Color Online) Competition among different sources for direct photon production in Au+Au collisions at 
^/sNN =200 GeV for different centralities. The leading order contribution (LO) from primordial elementary scatterings is 
plotted as dotted lines, thermal contribution as dash-dotted lines, fragmentation contribution as dashed lines and jet-photon 
conversion as solid lines. 
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Figure 10: (Color Online) The contribution to Raa from fragmentation and conversion, for different centralities. 



contribution from the QGP phase at p t > 1 GeV/c. In 
the current setting of hydrodynamic simulations at the 
RHIC energy, the mixed phase exists for a very long time 
(~ 8 fm/c). This contributes mostly at lower p t values. 
By combining the results shown in Figs. [9land[TTl con- 
tribution of thermal radiation from the QGP phase is 
dominant in the region 1 < p t < 4 GeV/c. This momen- 
tum window may provide information inside the hot and 



dense matter, e.g., the initial temperature at the center, 
which may not be reached directly by hadron spectr 
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Figure 11: (Color Online) (a) Thermal contribution to Raa from 0-10% (dashed line), 20-30% (solid line), and 40-50% (dotted 
line), (b) Fractions of thermal photon yields from the QGP phase, the mixed phase, and HG phase. 



VI. CONCLUSION 

We calculated the centrality-dependence of pt spectra 
for direct photons in Au+Au collisions at the RHIC en- 
ergy, based on a realistic data-constrained (3+l)-D hy- 
drodynamic description of the expanding hot and dense 
matter, a reasonable treatment of propagation of par- 
tons and their energy loss, and a systematic considera- 
tion of main sources of direct photons. In this study, 
four main sources are considered, namely, leading order 
(LO) contribution from primordial elementary scatter- 
ings, thermal radiation from the fluids, fragmentation 
from hard partons, and jet photon conversion (JPC). 
Similar work [4l| has been done before the appearence 
of the most recent data j38[. Our results agree nicely 
with the recent low p t data. 

The role of jet quenching in the high pt region of direct 
photons production has been checked via fragmentation 
photons and jet photon conversion sources. For these 
two sources, the suppression of the photon rate due to 
parton energy loss is significant in central Au+Au colli- 
sions, and becomes less important towards peripheral col- 
lisions, similar to the suppression for meson production. 
Since experimentally one may separate isolated photons 
(LO+JPC) and associate photons (fragmentation pho- 
tons), our prediction may be tested in the future. 

Considering the total yields of direct photons, the con- 
tribution from fragmentation and conversion are small, 
contributing between 5% and 10%. However, parton en- 
ergy loss plays nevertheless an important role: Without 
it, these second order effects would contribute signifi- 
cantly Without jet quenching, the nuclear modification 
factors Raa would depend visibly on the centrality of 
the collisions. A strong energy loss is actually necessary 
to get the centrality scaling of Raa in our calculation - 



a centrality scaling which has observed by the PHENIX 
collaboration. In this sense, properties of the bulk mat- 
ter affect the photon yields at intermediate values of pt, 
via the parton energy loss. 

The low pt region is totally dominated by thermal radi- 
ation, providing direct information about the bulk mat- 
ter. We find that Raa of photons at pt values below 
1 GcV/c is centrality independent. With increasing p t , 
the Raa for peripheral collisions drops much faster than 
the one for central scatterings. On the other hand, ther- 
mal photons are mainly emitted from the QGP phase at 
Pt >1 GeV/c even though the mixed phase and the HG 
phase occupy bigger space and longer time. So the differ- 
ent behavior of Raa for central and peripheral collisions, 
in the range 1 GeV/c < p t <3 GeV/c, manifests the fact 
that the plasma in central collisions is hotter compared 
to peripheral collisions. 

Still more investigation are needed for a precise char- 
acterization of the properties of the plasma via thermal 
photons. Besides, the elliptic flow of direct photons (es- 
pecially thermal photons) should provide more informa- 
tion of the plasma, which will be discussed elsewhere. 
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